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Abstract: We recently showed that a heavy quark moving sufficiently fast through a quark- 
gluon plasma may lose energy by Cherenkov-radiating mesons [1] . Here we review our previous 
holographic calculation of the energy loss in M = 4 Super Yang-Mills and extend it to 
Q\ longitudinal vector mesons and scalar mesons. We also discuss phenomenological implications 

for heavy-ion collision experiments. Although the Cherenkov energy loss is an 0(1/N C ) effect, 
a ballpark estimate yields a value of dE/dx for N c = 3 which is comparable to that of other 
mechanisms. 
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1. Introduction 

With the advent of the Large Hadron Collider (LHC) the field of heavy-ion collisions (HIC) 
enters a new era. The center-of-mass energy per nucleon in LHC collisions, ^/s NN ~ 5.5 TeV, 
is almost 30 times larger than that of the most energetic collisions at the Relativistic Heavy 
Ion Collider (RHIC). The highest temperature of the quark-gluon plasma (QGP) created in 
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RHIC experiments is approximately T RH ic — 2T C , with T c ~ 175 MeV the deconfinement 
temperature of Quantum Chromodynamics (QCD). Despite the large increase in the collision 
energy, this is expected to lead only to a moderate increase in the plasma temperature at 
the LHC [2], i.e. T LHC ~ (3 — 4)T C . 1 In contrast, high-energy partons originating from hard 
initial collisions will be copiously produced at the LHC. This will allow the study of quarks 
and gluons in the 100 GeV range, an order of magnitude larger than that at the RHIC. 

Experimentally, extremely valuable information is obtained by analyzing the energy loss 
of these energetic partons as they travel through the QGP. In order to use this information to 
learn about the plasma, a theoretical, quantitative understanding of the different mechanisms 
of parton energy loss is needed. Several such mechanisms have been previously studied, both 
in QCD itself [3] and in the context of the gauge/gravity duality [4]. 

A remarkable conclusion from the RHIC experiments [5] is that the QGP does not behave 
as a weakly coupled gas of quarks and gluons, but rather as a strongly coupled fluid [6]. 
Because of the moderate increase in the temperature and the logarithmic running of the 
QCD coupling constant, a qualitatively rather similar behaviour may be expected for the 
QGP at the LHC. This makes it particularly important to understand mechanisms of parton 
energy loss that may operate at strong coupling. We recently uncovered one such mechanism 
[1] whereby a sufficiently fast heavy quark traversing a strongly coupled plasma loses energy 
by Cherenkov-radiating in-medium mesons. 

The analysis in [1] showed that this mechanism takes place in all strongly coupled, large- 
N c gauge theory plasmas with a gravity dual. The argument is so simple that we reproduce 
it in section 2 for completeness. This section emphasizes the universality of the mechanism, 
since no reference to a specific model is necessary. 

Ref. [1] also performed a quantitative analysis in the simple example of a quark moving 
through the N = 4 super Yang-Mills (SYM) plasma. The quark Cherenkov-radiates both 
vector and scalar mesons. The rate of energy loss into the transverse modes of the vector 
mesons was calculated in [1], and again we reproduce it here for completeness. The vector 
mesons in question are massive, and thus they also possess a longitudinal mode. Here we 
extend the calculation of [1] and obtain the rate of energy loss into longitudinal vector mesons 
and scalar mesons. The result for the former is qualitatively similar to that for the transverse 
modes, whereas the result for scalar mesons displays some qualitative differences. 

Ref. [1] presented a rather preliminary exploration of the potential implications of these 
results for HIC experiments. Here we elaborate on that discussion and extend it to include 
possible implications of the new results presented in this paper. 

2. A universal mechanism of quark energy loss 

The reason that the mechanism we are going to describe is universal is that it only relies on 

: A rough estimate is obtained by assuming that the temperature scales as the fourth root of the energy 
density. 
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two universal features of the gauge/gravity duality: 2 (i) the fact that the deconfined phase 
of the gauge theory is described by a black hole geometry on the gravity side [7], and (ii) 
the fact that a finite number N { of quark flavours is described by N { D-brane probes [8] - see 
fig- 1. 




Figure 1: D-branes and an open string in a black hole geometry. 

In addition to the gauge theory directions, the gravity description always includes a radial 
direction (denoted by p in fig. 1) which is dual to the gauge theory energy scale. The radial 
position of the horizon is proportional to the plasma temperature T. The D-branes extend 
in the radial direction down to a minimum value proportional to the quark mass M q . 

As it is intuitively clear, for M q sufficiently larger than T the D-branes sit completely 
outside the horizon [9, 10, 11, 12, 13]. 3 In this phase, scalar and vector gauge theory mesons 
are described by small, normalizable fluctuations of scalar and vector fields propagating on 
the branes, whose low-energy dynamics is governed by a Maxwell-like theory. The spectrum 
of these fluctuations is discrete and gapped, which means that stable heavy meson states exist 
in the plasma. In other words, sufficiently heavy mesons survive deconfinement, in agreement 
with lattice and potential model predictions for QCD [14]. 

Consider now the in- medium dispersion relation uj(q) for these heavy mesons, where uj 
and q are the energy and the spatial three-momentum of the meson, respectively. As an 
illustrative example, the dispersion relations for vector and scalar mesons in the N = 4 SYM 
plasma 4 are depicted in fig. 2. 

As q — > oo, the DR becomes linear: u){q) ~ v hm q, with v hm < 1. This subluminal limiting 
velocity, which is the same for all mesons, can be easily understood in the gravitational 
description [10]. Since highly energetic mesons are strongly attracted by the gravitational 
pull of the black hole, their wave-function is very concentrated at the bottom of the branes. 
Consequently, their velocity is limited by the local speed of light at that point, v Xim (see 

2 In the limit N c ,gYy[N c — > oo. 

3 In contrast, as the ratio M q /T decreases, a first-order phase transition eventually occurs and a part of the 
branes falls through the horizon. See section 3. 

4 With quarks introduced as D7-branes; see below for details. 
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Figure 2: Left: Dispersion relation for the transverse (black, continuous curve) and longitudinal (red, 
dashed curve) n = modes of a heavy vector meson with w lim = 0.35 in the AT = 4 SYM plasma. 
Right: Analogous curves for a scalar (black, continuous curve) and pseudoscalar (red, dashed curve) 
meson. In both plots the blue, continuous straight lines correspond to u: = vq with w lim < v < 1. The 
black, dotted, vertical lines mark the crossing points between the meson dispersion relations and the 
blue lines. 



fig. 1). Because of the black hole redshift, v lim is lower than the speed of light at infinity (i.e. 
at the boundary), which is normalised to unity. In the gauge theory this translates into the 
statement that u Um is lower than the speed of light in the absence of a medium, namely in 
the vacuum. The reason is that the absence of a medium in the gauge theory corresponds to 
the absence of a black hole on the gravity side, in which case v Um = 1 everywhere. 

Imagine now a heavy quark in the plasma. In the gravitational picture, this is described 
by a string that starts on the D-branes and falls through the horizon - see fig. 1. In order to 
model a highly energetic quark we consider a string whose endpoint moves with an arbitrary 
velocity v at an arbitrary radial position p^. Roughly speaking, the interpretation of po 111 the 
gauge theory is that of the inverse size of the gluon cloud that dresses the quark. This can 
be seen, for example, by holographically computing the profile of (TrF 2 (x)} around a static 
quark source dual to a string whose endpoint sits at p = po [15]. 

Two simple observations now lead to the effect that we 

, . _,, „ • , , i ■ meson/ 

are interested in. Ihe first one is that the string endpoint 

is charged under the scalar and vector fields on the branes. quark cfuark 

In the gauge theory, this corresponds to an effective quark- 
meson coupling (see fig. 3) of order e ~ l/y/N~ c . We will Figure 3: Effective quark- 
derive these facts rigorously below, but physically they can meson coupling, 
be understood very simply. The fields on the branes describe 

fluctuations around the branes equilibrium configuration. The string endpoint pulls on the 
branes and therefore excites (i.e. it is charged under) these fields. The branes tension is of 
order l/g s ~ -/V c , where g s is the string coupling constant, whereas the string tension is N c - 
independent. This means that the deformation of the branes caused by the string is of order 
e 2 ~ l/iV c . We thus conclude that the dynamics of the 'branes+string endpoint' system 
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is (a generalization of) that of classical electrodynamics in a medium in the presence of a 
fast-moving charge. 

The second observation is that the velocity of the quark may exceed the limiting velocity 
of the mesons, since the redshift at the position of the string endpoint is smaller than at 
the bottom of the branes. As in ordinary electrodynamics, if this happens then the string 
endpoint loses energy by Cherenkov-radiating into the fields on the branes. 5 In the gauge 
theory, this translates into the quark losing energy by Cherenkov-radiating scalar and vector 
mesons. The rate of energy loss is set by the square of the coupling, and is therefore of order 
1/N e . 

3. Quarks in the J\f = 4 SYM plasma 

The four-dimensional N = 4 SYM theory with gauge group SU(N C ) at non-zero temperature 
is dual to type IIB string theory on the gravitational background sourced by N c black D3- 
branes. N f quark flavours may be introduced in the gauge theory by adding D7-branes on 
the gravity side. The relative orientation of the 'colour' and 'flavour' branes is summarised 
by the array 

0123456789 
N c D3: x x x x (3.1) 

N { D7: xxxxxxxx 

In the limit N { <C N c the backreaction of the D7-branes on the spacetime metric may be 
ignored and the D7-branes may be treated as probes in the gravitational background sourced 
by the D3-branes. Following [1, 10] we write the spacetime metric as ds 2 = L 2 ds 2 (G), where 



ds\G) = £ 



Cdt 2 + fdx 2 



+ [dr 2 + r 2 dn 2 3 + dR 2 + R 2 d$ 2 ] , (3.2) 



and 



L 4 = 4Trg s NJ 4 s , p 2 = R 2 + r 2 , / = 1 - 1/p 4 , f = 1 + 1/p 4 . (3.3) 

The four gauge theory directions are x M = {t, x} = {t,x 1 }, and they are identified with the 
0123-directions shared by both sets of branes in (3.1). The metric inside the second set of 
brackets in (3.2) is just the flat metric on R 6 = I 4 x I 2 , which corresponds to the 456789- 
directions in (3.1), written in terms of two sets of spherical coordinates {r, f^} and {R, 
The coordinate p is the overall radial coordinate in M 6 . This splitting is convenient since the 
D7-branes extend along the {r, f^j-directions. All coordinates above are dimensionless, and 
they are related to their dimensionful counterparts (denoted with tildes) through 

a? = irTx», {r,R,p} = -^ f {f,R,p}. (3.4) 

In particular, this means that we are measuring energy and momentum in the gauge theory 
in units of ttT. In addition, since the horizon of the metric (3.2) in dimensionless coordinates 



3 This can be viewed as a particular limit of string breaking - see sec. 7. 
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lies at p hm = 1, we see that the size of the horizon in physical units is proportional to the 
gauge theory temperature, i.e. p hor oc T. 

D7-brane embeddings 

We now specialize to N { = 1; we will discuss the case N f > 1 in sec. 7. The action gov- 
erning the dynamics of a D7-brane in the background sourced by D3-branes takes the form 



-T D7 J d 8 aV- det (g + 2tt£ 2 F) + T D7 J C 4 A F A F 



(3.5) 



In this equation T D7 = l/g s (2n) is the D7-brane tension, £ s is the string length, x a (a = 
0, . . . , 7) are intrinsic coordinates on the brane's worldvolume, g is the induced metric, F = dA 
is the field-strength of the worldvolume U(l) gauge field A a , and C4 stands for the pull-back 
of the spacetime Ramond-Ramond four-form potential sourced by the D3-branes. As we will 
see below, the term in the action involving C4 will not contribute to any of our calculations. 

In order to describe the D7-brane embedding we use x a = {x^, r, O3} as worldvolume co- 
ordinates. In other words, the brane extends along the gauge theory directions and the radial 
direction r, and it wraps an S 3 in the directions transverse to the D3-branes. Translational 
symmetry along and rotational symmetry along then imply that the embedding must 
be specified as R = R(r) and 1? = i?(r). Since 1? is also a symmetry direction, a consistent 
solution is obtained by choosing $ = const. A typical D7-brane embedding with different sets 
of coordinates suppressed is shown in fig. 4. 



M„ tx 




D7-brane 



M, x 




Figure 4: Minkowski-type D7-brane embedding showing the S 3 wrapped by the branes (left) and the 
gauge theory directions (right). The angular coordinate 1? is suppressed in both cases. The asymptotic 
position of the brane is proportional to the quark mass M q , whereas the size of the black hole horizon 
(shown in dark grey) is proportional to the temperature T . 

Under these circumstances the induced metric on the D7-brane takes the form ds 2 = 
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L 2 ds 2 (g) with 



ds\g) = 



f 2 jj.2 , 7 j->1 



f 



dt z + fdx z 



3 • 



(3.6) 



where i? = dR/dr. The function i?(r) is determined by inserting (3.6) in (3.5), setting F = 0, 
and varying with respect to R(r). The resulting Euler-Lagrange equation of motion is 



1 



1 



(r 2 + R 2 rJ 7lT^2 



R 



r A R 



(r 2 + i? 2 ) 5 



V 7 ! +« 2 



In the limit r — > oo, this equation leads to the asymptotic behaviour 

. c 
R(r) ~ m + -~ H . 



(3.7) 



(3.8) 



Holography relates the dimensionless constants m and c to the quark mass and condensate 
as (see [10] for details) 



M q = -V\Tm, {^) 



1 



\N { N c T 6 c, 



(3.9) 



where A = 5f 2 M A^ c = 2irg s N c is the 't Hooft coupling. An important point is that the constant 
m can also be written in terms of the mass M mcs of the lightest meson in the theory at zero 
temperature as [10]: 

2M q M mos 

m = ^T^ = ^F- ( 3 - 10 ) 



\/AT 2vrT 

Eq. (3.7) cannot be solved analytically, but numerical solutions for any value of the 
asymptotic brane position, m = R(r — > oo), can be easily found. The constants m and 
c correspond to the two solutions at infinity of the second-order equation of motion (3.7). 
These solutions are mathematically independent, but not physically: once m is specified, the 
requirement of regularity in the interior determines c. The physical solution is thus uniquely 
characterized by the value of m. 6 In the gauge theory this translates into the statement that 
once the quark mass (and the temperature) are specified, the dynamics determines the quark 
condensate. 

Solutions of eq. (3.7) fall into two classes. For m > 1.3, i.e. for quark masses sufficiently 
larger than the temperature, the brane bends towards the horizon because of its gravitational 
pull, but the brane tension is able to compensate for this and the brane sits entirely outside the 
horizon, as in fig. 4 and on the left-hand side of fig. 5. In this case the brane embedding is of 
the so-called 'Minkowski' type [9, 10], and we will denote by Rq = R{r = 0) the radial position 
of the bottom of the branes. For Minkowski embeddings there is a one-to-one correspondence 
between m and Rq, and the lowest value m = 1.3 corresponds to Rq = 1.2. 



For thermodynamically stable embeddings. In the case of thermodynamically metastable or unstable 
embeddings, c may be multivalued [9, 10]. 
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Minkowski embedding 



Black hole embedding 




Figure 5: First-order phase transition between Minkowski and black-hole type embeddings. 

In contrast, for m < 1.3, the brane falls through the black hole horizon in a so-called 
'black- hole embedding', since in this case the induced metric on the branes possesses a horizon 
- see fig. 5. At m = 1.3 a first-order phase transition between the two phases occurs, as 
indicated in fig. 5. As we will see below, in the Minkowski phase stable mesons exist, and their 
spectrum is discrete and gapped. In contrast, no stable mesons (in fact, no quasi-particles) 
exist in the black hole phase. For this reason, the phase transition above is interpreted in the 
gauge theory as a dissociation or 'melting' phase transition for mesons in the Af = 4 plasma 



The value m = 1.3 thus corresponds to the (stable) Minkowski- type brane that comes 
closest to the horizon, and therefore to the one for which the in-medium meson dispersion 
relation is most dramatically modified with respect to that in the vacuum. For this reason, 
we have chosen the embedding with m = 1.3 to illustrate some results in the sections below. 
Specifically, we see from eq. (3.6) that the local speed of light at the bottom of the brane is 



where we recall that Rq = R(r = 0). Since m = 1.3 corresponds to Rq = 1.2, the formula 
above gives u Um ~ 0.35, i.e. in this case the limiting velocity of mesons in the plasma is about 
1/3 of that in the vacuum. 

4. Meson dispersion relations in the N = 4 SYM plasma 

Despite the fact that the M = 4 SYM theory is in a deconfined phase at any T > 0, stable 
quark-antiquark states exist for sufficiently large M q /T, and the spectrum of these mesons 
is discrete and gapped. In particular, scalar and vector mesons in the gauge theory are dual 
to regular, normalizable modes of the scalar and vector fields on the D7-brane. Here we will 
review the dispersion relations for these modes, which we will need in order to compute the 
quark energy loss below. The dispersion relation for (some) vector mesons in the D3/D7 



[9, 10, 16]. 




(3.11) 
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system appeared in [1], but no details were presented there. The dispersion relation for scalar 
mesons was first computed in [10] and then revisited in [17]. Here we will review the result in 
the geometric parametrization of [17], which is particularly suited for calculating the energy 
radiated into these modes by the quark. 

4.1 Vector mesons 

These are dual to regular, normalizable fluctuations of the worldvolume gauge field A. The 
AT = 4 SYM theory possesses an internal, global SO (6) symmetry that is broken down to 
6*0(4) by the addition of quarks. In the string description, SO (6) is the isometry group of 
the spacetime metric (3.2), whereas SO (A) is the isometry group of the £ 3 wrapped by the 
D7-branes. Under the preserved £0(4) symmetry, meson modes decompose into singlet and 
non-singlet modes. Since we are interested in using the N = 4 SYM plasma as a toy model for 
the QCD plasma, and since QCD possesses no analog of the £0(4) symmetry, we will focus 
on singlet modes. The equation of motion for these modes receives no contribution from the 
second term in the action (3.5) [18], and therefore we will ignore this term in the following. 

In conclusion, since we are interested in singlet modes, we only need to consider the first 
term in the action (3.5). In addition, since we are only interested in their dispersion relation 
(as opposed to higher-order couplings), it suffices to expand this term to quadratic order in 
F in the fixed worldvolume metric (3.6). The result is 

£v CCtor = -T D7 L 4 (27r/ s 2 ) 2 J d s x^-g \F ab F ab , (4.1) 
leading to the equation of motion 

^ VaF ab = da{ ^ F ab )=0 ^ (42) 

The metric g that enters these expressions is that in eq. (3.6), which contains no factors of 
L; these have been explicitly included in the prefactor of (4.1). 
Singlet modes take the form 

A ll = Ap{xP,r), A r = A r (x»,r), An 3 = , (4.3) 

i.e. they have no components along the £ 3 and depend only on r and the gauge theory 
directions. The equations of motion are further simplified by the gauge choice A r = 0, which 
we will employ henceforth. In addition, we will work with the Fourier components of the 
gauge field defined through 

M*> *.»■)=/ M»> 9. *) e ~ m ' X . ( 4 - 4 ) 
J (2tt) 

where oj and q are the energy and the three-momentum of the meson, respectively. Finally, 
we choose the momentum to point along x 1 without loss of generality. 



7 At higher orders gauge field fluctuations would mix with scalar fluctuations. Similarly, at higher orders 
singlet modes would generically mix with non-singlet modes. 
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Under the conditions above, the equations of motion for the transverse modes A2,A^ 
decouple from each other and from those for the longitudinal mode Aq,Ai, so we will study 
them in turn. 



Transverse modes 

Let us collectively denote A = {A2, A3}. Both modes obey identical equations of motion 
which take the form 

d r {V9 9 rr g 33 d r A) - ^gg 33 (gV + g n q 2 ) A = 0. (4.5) 
Upon using (3.6) this becomes 

d r (-^5^) + 7TT^4 (^-^]a = 0. (4.6) 

Since we are interested in regular, normalizable solutions, we can expand A as 

A(u>, q,r)=J2 A»(w, q) £„(g, r) (4.7) 
n 

in terms of a basis of regular, normalizable eigenfunctions {Cn(q,r)} in the radial direction. 
These are solutions of eqn. (4.6) with g-dependent eigenvalues oj = uj n (q), he. they obey the 
eigenstate equation 

-d r [ d r UQ, r) J + Vl + R 2 i^- q 2 UQ, r) = Vl + R 2 ^ ^(qf Uq, (4-8) 

and satisfy the orthonormality relations 

/ dr J —Vl + R 2 U{q,r)Uq,r) = 6 mn . (4.9) 
Jo JP 

As we will see in more detail below, the discreteness of the spectrum is guaranteed by the 
boundary conditions on the £„: regularity at r = and normalizability at r = 00. Inserting 
the expansion (4.7) in (4.6), and using the eigenstate equation (4.8) and the orthonormality 
relations (4.9), we find that each of the A n (uj,q) fields obeys an independent equation of the 
form 

[u J 2 -u J 2 n (q)}A n (u J ,q)=0. (4.10) 

Thus, through the expansion (4.7) we have 'Kaluza-Klein-reduced' the five-dimensional field 
A(u, q, r) to a discrete, infinite tower of independent four-dimensional fields {A n (uj, q)}. Each 
of these fields is dual in the gauge theory to a transverse vector meson with dispersion relation 
oj = u n (q), which is the physical meaning of the wave equation (4.10). In the gauge theory, 
each of the mesons in this infinite set is distinguished by its 'internal' quantum number n. In 
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Figure 6: Normalized transverse vector meson radial profiles £,n{q,r) for a D7-brane embedding with 
to = 1.3. The blue, violet, brown and green curves (i.e. bottom to top on the left, top to bottom on 
the right) correspond to q — 1,2,4, 11, respectively. The rescalings on the right-hand side correspond 
to those in appendix C. Note that the area under the curves is not unity because of the non-trivial 
measure in eq. (4.9). 

the string description, each value of n corresponds to a different, (/-dependent radial 'wave- 
function' Cn(Q, r )- As we will see below, this structure of mesons in the fifth dimension will 
play an important role in determining the strength with which each of them couples to a 
quark. 
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Given that the brane embedding R(r) entering eq. (4.8) is only known numerically, the 
radial profiles must also be found numerically. The general solution of eq. (4.8) behaves as 
£ n ~ a + b/r 2 as r — > 0, and as £ n ~ a+b/r 2 as r — > oo, for some constants a, b, a, b. Regularity 
at r = requires 6 = 0, whereas normalizability imposes the condition a = 0. For fixed q, 
these two requirements are mutually compatible only for a discrete set of values of the energy, 
ui n (q). This is the origin of the dispersion relation. 

Fig. 6 shows several numerically-obtained 
radial profiles of the first lowest-lying modes 
£,n(Q, r ) for several values of q for a D7-brane 
embedding with m = 1.3. The corresponding 
values of the energy, w n (g), are given by the 
dispersion relation curves in figs. 2 and 7. As 
is familiar with solutions of Schrodinger-like 
equations, the nth solution possesses n zeros. 
More importantly, we see that the radial wave- 
functions for all these modes become concen- 
trated around the bottom of the brane, r ~ 0, 
as q — > oo. Relatedly, we observe that the lim- 
iting velocity of all these modes agrees with 
the local speed of light at the bottom of the 
brane, eq. (3.11), as expected from the general 
argument in section 2. 




Figure 7: Dispersion relation for the first trans- 
verse modes £„ of a heavy vector meson with 
vum = 0.35 in the TV = 4 SYM plasma. The 
curves correspond to n = 0,1,2 from bottom 
to top. The blue, continuous straight line cor- 
responds to lu = vq with i; lim < v < 1. 



Longitudinal modes 



Eq. (4.2) with b = 0, 1 and r yields two second-order dynamical equations and a first-order 
constraint equation, respectively, in which the longitudinal components Aq and A\ are cou- 
pled to one another. Only two out of the three equations are independent, which we take to 
be 



Or 



0. 

o. 



-g g "g^d r Ao)+iq^9g n 9 00 E 
icog 00 d T A - iqg n d r A 1 

where we have introduced the gauge-invariant electric field 

E = Fio = iqAo + iujA\ . 
From these two equations, it is easy to see that E satisfies 

-iqd r o 11 , — Y~m drE )+ l lV-9 9 9 E 
V 19 + w 9 J 

In order to turn this into an eigenstate equation we introduce a new field $ defined as 



— rr 00 11 

-99 9 9 



q 2 g 11 + w z g' 



2 n 00 Wr 



d r E. 



(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 
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Eq. (4.14) then implies the inverse relation 

E= -== \ 1 00 a r $ (4.16) 

which, when substituted back into (4.15), yields the equation of motion for $: 

/ 1 \ q 2 q U + w 2 q m 

Inserting the explicit form of the metric functions, we arrive at 

dr ( . . dA + 2^^ p - £^ * = . (4.18) 



From this point onward, we proceed as in the case of transverse modes. We expand as 

g, r) = $ n (w, g)0„(g, r) (4.19) 

n 

in terms of a basis of regular, normalizable eigenfunctions {0n(<7j r)} in the radial direction. 
These are solutions of eq. (4.18) with q-dependent eigenvalues u = uj n (q), and are subject to 
the orthonormality relations 



poo 



/ dr2 —3 — cf) n (q,r)<p m (q,r) = 6 mn . (4.20) 

Jo J r 

As in the case of transverse modes, the longitudinal modes 3> n (w, q) obey the wave equation 

[uj 2 - co 2 n (q)} $„(0J,q) =0, (4.21) 

as appropriate for a four-dimensional field with dispersion relation oj = uo n (q). Again, through 
the expansion (4.19) we have Kaluza-Klein-reduced the five-dimensional field $>(u),q,r) to a 
discrete, infinite tower of independent four-dimensional fields {ip n (oj,q)}, each of which is 
dual to a longitudinal vector meson in the gauge theory. 

The general solution of eq. (4.18) behaves as <j) n ~ a + br 4 as r — > 0, and as 4> n ~ a+61ogr 
as r — > co, for some constants a, b, a, b. Normalizability with respect to (4.20) as r — > requires 
that a = 0, and regularity as r — > 00 requires that 6 = 0. As in the case of the transverse 
modes, for fixed q these two requirements are mutually compatible only for a discrete set of 
energies u n {q). 

Fig. 8 shows the electric field, eq. (4.16), of several numerically-obtained radial profiles 
of the lowest-lying mode (f) n =o(q,r) for several values of q. The corresponding values of the 
energy, w n= o(q), are given by the dispersion relation curve in fig. 2. These results correspond 
again to a D7-brane embedding with asymptotic position m = 1.3. We observe the same 
limiting velocity t> lim = 0.35 given by the local speed of light at the bottom of the branes, 
eq. (3.11). 
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Figure 8: Electric field E (q,r) of the normalized longitudinal vector meson radial profiles (j>o(q,r) 
for a D7-brane embedding with m = 1.3. The blue, violet, brown and green curves (i.e. bottom to 
top on the left, top to bottom on the right) correspond to q — 1,2,4, 11, respectively. The rescalings 
on the right-hand side correspond to those in appendix C. Note that the area under the curves is not 
unity because of the non-trivial measure in eq. (4.20). 

4.2 Scalar mesons 

The scalar fields on the brane get excited by the string endpoint because the string tension 
pulls on the brane. A crucial feature is the fact that the boundary conditions at the string 
endpoint imply that the string must end orthogonally on the brane. (The unfamiliar reader 
can find a concise derivation in appendix A.) For this reason it is convenient to work with 
spacetime coordinates that locally parametrize the directions orthogonal to the brane. In the 
case of interest to us, the coordinate in (3.2) satisfies this requirement, since the vector field 
d/d"d is orthogonal to the fiducial D7-brane embedding at each point on the brane. However, 
the R coordinate does not meet this requirement, since d/dR is in general not orthogonal 
to the brane due to the brane bending in the r — R directions. We therefore follow [17] and 
work with two geometric coordinates X A defined as follows. At each point on the brane, the 
two-dimensional space orthogonal to the brane is spanned by the unit vectors 

55 -*■•>!■ v ^Tr (4 ' 22) 

so a general vector orthogonal to the brane takes the form U = X a Va- For each vector U we 
shoot off a geodesic with unit affine parameter that at the brane has U as its tangent vector. 
The endpoint of this geodesic defines a point in a neighborhood of the brane. In this way we 
obtain a one-to-one correspondence (the so-called exponential map [19]) between the values 
of X A and the points near the brane. In other words, X A are the coordinates orthogonal to 
the brane that we were seeking, since on the brane we have d/dX A = Va by construction. In 
particular, this implies that in these coordinates Gab = &AB when evaluated precisely on the 
brane, where G is the spacetime metric (3.2). Note that in this section the ten dimensionless 
coordinates of spacetime are thus {x a , X A }, with a = 0, ... ,7 and A = 1,2. 

We chose the fiducial embedding of the brane to be given by X A = 0, so that the X A 
fields parametrize fluctuations around it. As shown in [17], to quadratic order in these fields 
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the D7-brane action takes the simple form 



^scalar 



7 



T D7 L 8 / cPxy/^g 



l -g ab d a X A d b X B G AB + l -m 2 AB (x)X A X B 



(4.23) 



where g is the induced metric (3.6) on the fiducial embedding of the brane, i.e. it is X A - 
independent. As usual, this metric contains no factors of L, since this have been factored out 
explicitly in front of the action. The position-dependent mass matrix m 2 AB {x) is diagonal and 
given in terms of geometric quantities as 

m 2 n = Rn + #2112 + 2# 22 + - R , 

m 2 22 = -R 22 + #2112 , (4.24) 

where 

#2112 = V^VfVfV^RMNPQ , (4.25) 

#n = V^VfRMN , (4.26) 

#22 = V 2 M V 2 N R MN • (4.27) 

Rmnpq and Rmn are, respectively, the Riemann and the Ricci tensors of the ten-dimensional 
spacetime metric G, R is the corresponding Ricci scalar, and ( 8 )# is the Ricci scalar of the 
eight-dimensional induced metric on the brane g. Again none of these quantities contains any 
factors of L. 

Using the fact that Gab = &AB the action (4.23) leads to the equation of motion 

JZTg V 2 X A - m 2 X A = d a ( ^~g d a X A ) - ^~g m 2 X A = 0, (4.28) 

where m = m\\ or 771,22 as appropriate. Following the vector meson case, we focus on the 
zero-mode of X A on the S 3 , and work with its Fourier components X a (lj, q, r), for which the 
equation of motion is 

d r (^~9 9 rr d r X A ) - ^g ( 5 °V + g n q 2 + to 2 ) X A = , (4.29) 

which upon substitution of the metric functions becomes 

d r (^£LdrX A ) + f]r 3 ^¥ (%u 2 - *f — to 2 ^ X A = . (4.30) 

Wi + # 2 / \pf p 2 f J 

As usual, we expand X A as 

X A (u, q,r)=Y J X A (u, k) ip A (k, r) (4.31) 

n 

in terms of a basis of normalizable eigenfunctions {ip A (k,r)} in the radial direction. These 
are solutions of eq. (4.30) with (/-dependent eigenvalues ui = LO A (k), and are subject to the 
orthonormality relations 

/ dr 2/ ^i(q,r)^ A (q,r) = 6 mn . (4.32) 

Jo PI 
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Inserting the expansion (4.31) in (4.30), and using the orthonormality relations (4.32), we 
find that each of the X a (lj, (?) fields obeys 

[«, 2 -ul A (q)]x£(u>,q)=0, (4.33) 

as expected. As explained in [12], the modes X a (uj, q) with A = 1,2 correspond in the gauge 
theory to scalar and pseite/oscalar mesons, respectively. 

Both masses (for A = 1,2) in eq. (4.30) behave as m 2 ~ —3 — m 2 /r 2 + • • • for r — > oo 
and m 2 ~ —c\ + C2r 2 + • • • for r — > 0, where c\ t 2 are positive constants. It follows that the 
two independent solutions for X A behave as 1/r and 1/r 3 for r — > oo and as r° and 1/r 2 for 
r — > 0. Thus in this case normalizability requires that X A ~ 1/r 3 for r — >■ oo and regularity 
requires that X" 4 ~ r° for r — > 0. As in the case of vector modes, for fixed q these two 
requirements are compatible with each other only for a discrete set of energies u A (q). 

Fig. 9 shows several numerically-obtained radial profiles of the lowest-lying mode <p n= o(q, r) 
for several values of q. The corresponding values of the energy, oj n= o(q), are given by the dis- 
persion relation curves in fig. 2. These results correspond again to a D7-brane embedding 
with asymptotic position m = 1.3. We observe the same limiting velocity u lim = 0.35 given 
by the local speed of light at the bottom of the brane, eq. (3.11). 

5. Quark energy loss in the Af = 4 SYM plasma 

We now turn to the main topic of this paper, namely the rate at which a heavy quark 
traversing the N = 4 SYM plasma loses energy by Cherenkov-radiating mesons. As we will 
see in detail below, the quark acts a source for the brane stress-energy tensor, defined as 

2 6S m 

This leads to its non-conservation, V a T a b / 0, where V is the covariant derivative defined by 
the eight-dimensional worldvolume metric g. Under these circumstances the energy per unit 
time deposited on the brane by the quark is given by 

— = -J d 7 xV=9~V a T a0 , (5.2) 

where the integral is taken over the brane's worldspace. For the reader's convenience, a short 
derivation of this formula is provided in appendix B. Our task below will be to evaluate this 
formula for the cases of vector and scalar mesons. 

5.1 Energy loss into vector mesons 

The endpoint of an open string attached to the brane couples to the worldvolume gauge field, 
so the action (4.1) is modified in the presence of the quark to 

S vcctoI = -T D7 L 4 (27Tl 2 Y J d 8 x^g- -F ab F ab - J drA a — , (5.3) 
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Figure 9: Normalized scalar (top plots) and pseudoscalar (bottom plots) meson radial profiles <p n (q, f) 
for a D7-brane embedding with m = 1.3. The blue, violet, brown and green curves (i.e. bottom to 
top on the left, top to bottom on the right) correspond to q — 1,2,4, 11, respectively. The rescalings 
on the right-hand side correspond to those in appendix C. Note that the area under the curves is not 
unity because of the non-trivial measure in eq. (4.32). 



where, as usual, the worldvolume metric g contains no factors of L. The second term is 
the familiar coupling between an Abelian gauge field and a charged particle moving along a 
worldline parametrized as x a (r). 8 In order to work with a canonically normalized gauge field 
we rescale A — > eA with 

1 8^ 4 



e 2 = — = , (5-4) 

T D7 (2irP s ) 2 L 4 N c 



so that the action becomes 



S vcctor = - J d 8 x^g-F ab F ab - , j ,r.rA„.r 

where J a = 5^ 7 '(x — x(r)) x a . As anticipated in section 2, the coupling e between the quark 
and the vector mesons is of order 1/yiVc, which justifies our neglect of terms of order higher 
than quadratic in the action. 



J d s xA a J a , (5.5) 



8 The relative normalization between the two terms in the action can be confirmed by noting that it ensures 
that supersymmetric Blon-like excitations on a D7-brane in flat space have tension 1 /2-k£ 2 s , as in [20] . 
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Using the definition (5.1), the contribution from the gauge field to the brane's stress- 
energy tensor is easily calculated to be 

T ab = F ac F b c - l A g ab F 2 . (5.6) 

In the presence of the string endpoint, the equation of motion (4.2) for the gauge field is 
modified to 

^~gV a F ab = d a {^-gF ab ) = ej\ (5.7) 
which implies the non-conservation of the stress-energy tensor 

V=9~V a T ab = eF ba J a . (5.8) 

Inserting this into the general formula (5.2) yields the rate at which the quark deposits energy 
into the gauge field: 

= ~ej d 3 xdrdn 3 F 0a J a . (5.9) 



dE v 



dt 

This formula has the simple interpretation of minus the work done on the quark by the gauge 
field. In order to evaluate it, we need to specify the quark trajectory. For simplicity, we will 
assume that the quark moves with constant velocity along a straight line at constant radial 
and angular positions, so we write 

J a = 5 {3) {x-vt)5(r-r )5 {3) {n-n ) x (1,^,0,0). (5.10) 

In reality, ro and v will of course decrease with time because of the black hole gravitational 
pull and the energy loss. However, we will concentrate on the initial part of the trajectory 
(which is long provided the initial quark energy is large) for which ro and v are approximately 
constant [21] - see fig. 1. The delta- functions in (5.10) allow us to perform the integral in 
(5.9) and obtain 

dF, 

— ^ = -ev l F 0l (t, vt, r , Q ) . (5.11) 

We thus see that we need to compute the electric field sourced by the string endpoint at 
the location of the string endpoint itself. To do so, we will solve the equation of motion (5.7) 
by expanding the gauge field in normalizable modes in the radial direction, as in section 4. 
Note that the fact that the quark is localized on the S 3 means that it will radiate both into 
S* 3 singlets and non-singlets. A simple group theory argument shows that these two types of 
contributions can be calculated separately and independently at the quadratic level at which 
we are working. For the reasons explained in section 4, we will only calculate the energy loss 
into singlet modes, whose form in Fourier-space we recall to be: 

A^ = A^(u,q,r), A r = , Aq 3 = . (5.12) 

Without loss of generality, we choose q = (q, 0, 0) and v = (v cos 0, v sin 9, 0). After integrating 
over the S 3 , the relevant Fourier-space components of the current are then 

= 2ir5(u-qv cos 6) 5(r - r ) x (l,i;cos0,i;sin0,O). (5.13) 
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We are now ready to compute the energy loss into the transverse and longitudinal modes of 
the gauge field. 



Transverse modes 



With the choice above the only transverse mode of the gauge field excited by the source 
is A = A2, which couples to J = J 2 . The equation of motion (4.6) for this mode now 
becomes 

B (jp±) + V^A^(^-ff]A^J, (5.4, 

where e = e/ils and the volume factor = 2n 2 comes from integration over the 5 3 . We now 
follow section 4 and solve (5.14) by expanding A as in (4.7), where the radial eigenfunctions 
{£n(<2S r )} satisfy exactly the same properties as in that section. In this case, inserting the 
expansion (4.7) in (5.14), and using the eigenstate equation and the orthonormality relations, 
we find that eq. (4.10) becomes 

[u? - c4(<?)] A n (u, q) = ej n (u, q) , (5.15) 

where 

J n (u,q)= / drj(u,q)£ n (q,r) = 2tt S (uj - qv cos9) v sind £ n (q, r ) . (5.16) 

An important fact implied by eqs. (5.15)-(5.16) is that each of the four-dimensional meson 
modes A n (oj,q) couples to the quark with an effective strength proportional to the value of 
£ ra at the location of the quark: 

e eS (q,ro) = e £„(<?, r ) . (5.17) 

The intuition behind this is that the radial profiles £n(<7,r) roughly play the role of a 'wave 
function' in the fifth dimension for the corresponding meson mode A n {oJ,q)- This fact will 
play an important role below. 

With retarded boundary conditions, as appropriate for the reaction to the quark's pas- 
sage, the solution of eq. (5.15) is 



An{u, q) = —2 • (5.18) 



ej n (uj,q) 
{uj + ief - ul(q) 

In order to evaluate the energy loss (5.11), we first express i*02(i, vt, ro, fio) as an integral over 
its Fourier components: 9 



"-E'trans 9 



dt 



ev d t A 2 (t,vt,ro,tto) 
dud^q 



Note that the singlet mode is independent of the S position f2o of the string endpoint. 
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Inserting the expansion (4.7) we obtain 
dE trans sr^ f dud 3 



dt 



E 



(2*)' 



ev sin 9){zu) An(u), q)£ n {q, r Q ) e 



-iojt gitq-v 



(5.20) 



Substituting the solution (5.18) for A n and using the delta- function in (5.16) to integrate over 
frequencies, we arrive at 



dF, 

trans 

dt 



E 



d 3 q 

W) 



{ev sm 9) {iqv cos 8) — — — - 2 — £ n {q, r ) 

{qvcosa + ie) — u^{q) 



(5.21) 



Note that the two exponentials have cancelled out upon setting cj = qvcosO. In order to 
integrate over momenta we set d 3 q = 2irq 2 dqdz, where z = cos 6, so that the integral above 
becomes 



dE tT . 



dt 



E 



e 2 v 



dq 
2^ 



dz 



z(l-z 2 ) 



_X 2tti (z + it) 2 - z 2 (q) 



(5.22) 




where z n (q) = v n (q)/v, and v n (q) = oj n {q)/q is the phase velocity of the re-th mode. 

The integral over z can be performed in 
the complex plane by considering the con- 
tour shown in fig. 10. The integral of interest 
corresponds to the integral over the segment 
T\. In the limit e — > 0, this coincides with 
the integral over T2, and the contribution 
from the vertical sides of the contour van- 
ishes. Thus in this limit the integral over 
z in (5.22) equals 1/2 times the contour in- 
tegral of fig. 10. Now consider the poles of 
the integrand, which lie at z = ±z n — ie. If 
v < v n (q) then z n (q) > 1 and the poles lie 

outside the contour, so the integral vanishes. In contrast, if v > v n (q), then both poles lie 
inside the contour and they yield identical contributions equal to (1 — z 2 ). Taking into con- 
sideration the extra minus sign coming from the orientation of the contour, the final result is 
thus 



Figure 10: Contour in the complex plane used to 
evaluate the z-integral in eq. (5.22). 



dE tI 



E 



e 2 v 



dt ^ 20 3 l n 2tt 

n 



dq 



<{q) 



1 



(5.23) 



We see that the energy loss is a discrete sum over all mesons, as well as an integral over all 
the momentum modes of each meson into which the quark is allowed to radiate. As expected 
for Cherenkov radiation, this can only happen if the velocity v of the quark exceeds the phase 
velocity of the corresponding momentum mode, v n (q). For example, in the case of transverse 
vector mesons, the quark can only emit momentum modes to the right of the dashed, vertical 
line in fig. 2. This cut-off is implemented by the Heaviside function in eq. (5.23). 
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Since the radial profiles £ n {q, fo) and the dispersion relations v n (q) entering eq. (5.23) are 
only known numerically, the energy loss must also be evaluated numerically. The result for 
the n = term in the sum is plotted in fig. 11. As one may expect, for fixed ro the energy 
loss increases monotonically with v up to the maximum allowed value of v, the local speed 
of light at tq. In other words, a quark sitting at a fixed radial position radiates more the 
higher its velocity is. As ro decreases, the limiting velocity of the quark approaches that of 
the mesons from above. Therefore the quark and the meson dispersion relation curves cross 
at a higher momentum, i.e. the vertical dashed line in fig. 2 moves to the right. This means 
that the characteristic momentum q chai of the modes contributing to the integral in (5.23) 
increases. As ro — > these modes become increasingly peaked at small r (see fig. 6(left)), 
and their effective couplings to the quark e eff (g char , ro) diverge. This explains why the energy 
loss at the maximum allowed value of the velocity diverges as ro — > 0. As we will discuss in 
section 7, however, this mathematical divergence is removed by physical effects that we have 
not taken into account. 




0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.88 0.90 0.92 0.94 0.96 0.98 1.00 
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Figure 11: Energy loss into the transverse vector mode Ar n —o) for an embedding with m — 1.32, R = 
1.2 (left) and m = 2.0, Rq — 2.0 (right). The continuous curves correspond (from top to bottom) to 
r = 0.86, 0.97, 1.10, 1.25, 1.45 (left) and to r = 1.50, 1.69, 1.91, 2.18, 2.52 (right). The dotted curve is 
defined by the endpoints of the constant-ro curves. 



Longitudinal modes 

In the presence of the source eq. (4.11) becomes 

d r {V=99 rr g 00 drA ) + iqy/^g n g 00 E = eJ° . (5.24) 
Eq. (4.12) remains unchanged, and together with the eq. above it yields 

d r f ^ri¥r44) d r E ) + i( i^~9 g ll g Q0 E = gj° . (5.25) 

V Q 9 + w g uu / 

We now introduce a new field $ defined as in eq. (4.15), but in this case the inverse relation 
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(4.16) is modified by the source: 

E = \ 1 0Q d r $ + . -=± n 0Q gJ° . (5.26) 

y / ZI gg 11 g w iqV^gg ll g m 

Substituting back into eq. (4.15) we find the new equation of motion for $: 



( ^koo **) + ^^ool * = br i^=^ ^ ) . (.1.27) 



Inserting the explicit form of the metric functions, we finally arrive at 
d r ( ? a,*) + 2^+^ U - CA * = -U( *f eA . (5.28) 



We now follow section 4 and solve (5.28) by expanding <3? as in (4.19), where the radial 
eigenfunctions {(f) n (q, r)} satisfy exactly the same properties as in that section. In this case, 
inserting the expansion (4.19) in (5.28), and using the eigenstate equation and the orthonor- 
mality relations, we find that eq. (4.21) becomes 

[u 2 - u 2 n (q)] <S> n (u, q) = eJ°(w, q) , (5.29) 

whose solution with retarded boundary conditions is 



eJn(u,q) 
(u + ief-uHq) 



The coefficients are given (after integration by parts) by 

/*°° 1 f 2tt 

Jn(u,q)= dr — - . J°(uj,q)d r (j) n (q,r) = —5(uj - qv cos 9)F n (q, r ) , (5.31) 

Jo ^/rVl + i? 2 n 

with 

Fn{q,r) = - 1 n Q Qd r (f) n (q,r) = ^ P . d r (p n (q,r) . (5.32) 

V 9 9 9 /r 3 v 1 + i? 2 

Note that the coefficients T n appear in the expansion of the electric field (5.26), i.e. 

E(u, q, r) = - V $ n (w, g)J" n (g, r) + - — ^= n QQ eJ°(w, g, r) . (5.33) 

iq v~9 9 9 

As in the case of transverse modes, we see from eq. (5.29) that the effective coupling between 
a longitudinal meson $ n (ui, q) and the quark is determined by the radial wave function of the 
meson, in this case 

e c a(q,r ) = eJ 7 n (g,r ) . (5.34) 

Our task now is to compute the rate of energy loss into longitudinal meson modes. For 
this purpose, eq. (5.11) instructs us to evaluate the electric field at the location of the quark. 
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If we naively do so using the expression (5.26) for E then the second term gives a divergent 
result, since J°(ro) oc 5(0). However, this divergence is unphysical: if one replaces the delta- 
function by a smooth charge distribution, then the integral over space in (5.9) vanishes. 
Indeed, suppose that the current (5.13) is replaced by 

J a = Q ( 3 ) (x - vt, r, n) x (1, v, 0, 0) , (5.35) 

where g is a smooth function. Then in Fourier space 

J a = 2tt5{uj - q-v) g i3) (q,r,n) x (1,^,0,0) (5.36) 

and the energy loss (5.9) is 

^ = — e J drdtt J d 3 x E(t, x,r,Q) ■ J(t, x,r,Q) 

= -ej drdU J ^e-^-^J ^ z E(u>, -q, r, fi) q, r, fi) , (5.37) 

where as usual z = cos is the relative angle between q and v and d 3 q = 2irq 2 dqdz. If we 
substitute the term in E that is proportional to J° we see that the integrand is proportional 

t0 1 

— S(u + qvz) 5(uj — qvz) g(—q, r, Q) g(q, r,fl) z . (5.38) 
iq 

This is odd under q — > —q (since u and cu are dummy variables) and therefore the integral 
over z vanishes. 

We therefore conclude that we can neglect the second term in (5.33) in order to evaluate 
(5.11). Following the previous section we have 



dEl °^ -= ev 1 E{t,vt,r ,tt Q ) 



dt 

Ef dojd?q 
, ~J w 



(ev cos 6) $ n (w, q)T n (q, r ) e^V*^ 



Ef d 3 q e 1 9/ \ 

— / tt, (evcosv) k J- n (q,ro) 

e 2 f°° dq 1 o . . Z* 1 dz z . 

^ y 2vr g i_! 2tu (z + ze) 2 - 4(g) 

The z-integral can be evaluated using the same contour of fig. 10. In this case each pole 
contributes -1/2, so the final result is 

^ = E^ -^0)0(1-^) . (5.40) 

dt ^ 2Q 3 v J 2ir q nW ' \ v J K J 

The energy loss into the n = mode is shown in fig. 12. The same comments as in the 
case of transverse modes apply here. 
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Figure 12: Energy loss into the longitudinal vector mode $( n =o) f° r an embedding with m — 
1.32, Rq — 1.2 (left) and m = 2.0, Rq = 2.0 (right). The continuous curves correspond (from top 
to bottom) to r = 0.86,0.97,1.10,1.25,1.45 (left) and to r = 1.50,1.69,1.91,2.18,2.52 (right). The 
dotted curve is defined by the endpoints of the constant-ro curves. 



5.2 Energy loss into scalar mesons 

The endpoint of an open string attached to the brane couples to the worldvolume scalar fields, 
so in the presence of the quark the action (4.23) is modified to 



5s, 



-T m L 
L 2 



<Fxy/=g 



l -g ab d a X A d b X B G AB + l -m 2 AB {x)X A X B 



j dry^GAB X A n B , 



(5.41) 



where x 2 = g ab x a x b and n A is the unit vector that is tangent to the string and orthogonal 
to the brane at the string endpoint. As usual, in the last term we have explicitly factored 
out the L-dependence associated to the spacetime metric (3.2). The square-root factor in the 
last term is necessary to make the integrand a scalar density on the worldline of the string 
endpoint, which is contained in the brane's worldvolume. We emphasize that the indices of 
x a are not contracted with the gauge theory metric rjab but with the D7-brane metric g a b- 
The relative normalization between the two lines in (5.41) can be confirmed as described 
in footnote 8. As explained in appendix A, the boundary condition at the string endpoint 
implies that the string ends orthogonally on the brane. The unit normal in the last term of 
the action (5.41) means that the string couples to the scalar that parametrizes the direction 
along which the string pulls on the brane, as one may intuitively expect. 

In order to work with canonically normalized fields we rescale X A — > \/T D7 L s X A so that 
the action becomes 



d 8 x, 



-9 



l -g ab d a X A d b X B G AB + l -m 2 AB (x)X A X B 



e / d 8 x J A X A , (5.42) 



where J a = 5^ (x — x{r)) \J — x 2 Gab n B and e is the same coupling constant defined in (5.4). 
Since the scalars do not interact with each other, they give independent contributions to the 
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stress tensor of the brane. Using the definition (5.1), the contribution from either scalar field 
is easily found to be 



1 



5m 



2 



r=-g V a T ab = JV b X- X 2 V b m + V a ( X 2 — b ) . (5.45) 



T ab = V a XV b X - -g ab [(VI) Z + m l X l \ + X 2 ^ . (5.43) 

For ease of notation, in this equation we have dropped the superindex 'A' on the scalar, and 
we will continue to do so below. The last term originates from the non-trivial dependence 
of the scalar masses (4.24) on the metric. Fortunately, we will see that we do not need to 
evaluate this term explicitly in order to compute the divergence of the stress-tensor. 

In the presence of the string endpoint, the equation of motion (4.28) for the scalars is 
modified to 

^~g (V 2 - to 2 ) X = e J . (5.44) 
Using this, the divergence of the stress tensor takes the form 

5m 2 ' 
5~g~ a 

The second term on the right-hand side is due to the possible spacetime dependence of the 
scalar masses m(x), but it vanishes identically for the case of interest here, 6 = 0, because of 
the time-translation invariance of the theory. The last term on the right-hand side vanishes 
when evaluated on a solution of the equations of motion, even in the presence of the source 
J. To see this, recall that the stress tensor must be identically conserved in the absence of 
the source because of the diffeomorphism invariance of the brane's worldvolume theory. This 
means, in particular, that when J = we have 

The key point now is that this equation is linear in X. Since the solution in the presence 
of the source is a linear supersposition of solutions of the source-less equation, linearity of 
(5.46) implies that this expression also vanishes for solutions of eq. (5.44) with J ^ 0. We 
thus conclude that (the time component of) the non-conservation of the stress-tensor in the 
scalar sector takes the form 

V=gV a T a0 = JV X . (5.47) 

As in the case of vector mesons, we consider a rectilinear quark motion with constant 
velocity, in which case 

J = ^-x 2 (r ) S^(x-vt) 5{r - r ) S (3) (0 - fi ) ■ (5-48) 

Note that, although the velocity v is the quark velocity as seen by a gauge theory observer, 
the prefactor above is not just a/1 — v 2 but depends non-trivially on the quark position in 
the radial direction through 

y/-x 2 (r ) = ^/-goo(ro) - gn(r )v 2 . (5.49) 
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Because of the black hole redshift, for fixed ro this factor vanishes before v reaches the speed 
of light, i.e. at v < 1. As we will see, this fact is responsible for a qualitative difference 
between the energy radiated into scalar and into vector mesons. 

Following section 4.2, we focus on the zero-mode of X on the S* 3 , and work with its 
Fourier components X(ui,q,r), for which the equation of motion takes the form 

d r (V=gg rr d r X) - (g 00 ^ + g n q 2 + m 2 ) X = eJ , (5.50) 

where the relevant Fourier-space components of the source are 

J = \J — x 2 {tq) 2tt5(uj — qv cos 6) 5(r — ro) , (5.51) 

and as usual e = e/Si^. We solve (5.50) by expanding X as in eq. (4.30). In the presence of 
the source, the equation obeyed by the X n mode is 

[u 2 - ul(q)] X n (u, q) = eJ n (u, q) , (5.52) 

where 

J n (uj,q) = J drJ(u,q)ip n (q,r) = A/-x 2 (r ) 2ir5(u) - qv cos 9)ip n (q, ro) ■ (5.53) 



X n (u,q)= , n r\,, - (5-54) 



With retarded boundary conditions, as appropriate for the reaction to the quark's passage, 
eq. (5.52) yields 

eJ n (uj,q) 
(w + ief - ul(q) 

As in the case of vector modes, we see from eq. (5.52) that the effective coupling between a 
scalar meson X n (oj, q) and the quark is determined by the radial wave function of the meson, 
in this case 

e eff (g, r ) = e ^/-x 2 (r ) ip n (q, r ) . (5.55) 

We are now ready to compute the rate of energy deposition into scalar mesons. From 
eqs. (5.2) and (5.47) we have 

rfffscalar = _ & J fafo^ £ J = _ & ^3^7^ ± (t ,Vt, V , Q ). (5.56) 

Following the steps of the vector meson case we find: 



Y^-ey/^W) J d ^{-iu)X n {u,q) Vn {q^)e- i ^ t ^ 

^2-e[-x 2 (r )] [ Jr^i-iqv cos 9)- 6 — <£(q,r ) 

„ J ( 27 ( qv cos 9 + is) -Lo 2 (q) 



(qv cos 9 + ie) — uj 2 (q) 
1 dz z 



= E^H 2 Mf g«rfto.n>)e(i-!^>) . (5 . 57) 
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The result for the energy loss into the lowest-lying scalar and pseudoscalar modes is shown 
in figs. 13 and 14. The main difference with respect to the case of vector mesons is the fact 
that the constant-ro curves do not rise monotonically as v increases, but instead they vanish 
when v reaches the local speed of light at tq. The reason for this is of course the factor in 
eq. (5.49). 
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V V 

Figure 13: Energy loss into the scalar mode ip^tz^ for an embedding with m = 1.32, Rq = 1.2 
(left) and m = 2.0, Ro — 2.0 (right). The continuous curves correspond (from top to bottom) to 
r = 0.86, 0.97, 1.10, 1.25, 1.45 (left) and to r = 1.50, 1.69, 1.91, 2.18, 2.52 (right). The dotted curve is 
defined by the maxima of the constant-ro curves. 
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Figure 14: Energy loss into the pseudoscalar mode <p) n=0 \ for an embedding with m = 1.32, Rq = 1.2 
(left-hand side plots) and m — 2.0, Rq — 2.0 (right-hand side plots). The continuous curves correspond 
(from top to bottom) to r = 0.86,0.97,1.10,1.25,1.45 (left) and to r = 1.50,1.69,1.91,2.18,2.52 
(right). The dotted curve is defined by the maxima of the constant-ro curves. 



6. Phenomenological implications for HIC experiments 

As is clear from our general discussion in section 2, the mechanism of Cherenkov energy loss 
depends only on two qualitative properties encoded in the dispersion relations of fig. 2: the 
fact that heavy mesons remain bound in the gauge theory plasma, and the fact that their 
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limiting velocity in the plasma is subluminal. Both properties can be motivated in QCD 
irrespectively of whether or not a string dual of QCD exists. The first property is suggested 
by the fact that sufficiently heavy mesons are smaller than the screening length in the plasma 
[22], and is supported by calculations of both the static quark-antiquark potential [23] and of 
Minkowski-space spectral functions in lattice-regularized QCD [24]. 10 The second property, 
which goes back to ref. [26], is suggested by the fact that moving mesons see a boosted, higher 
energy density that will melt them if they move sufficiently fast [27]. 11 

Rigorously verifying these two properties in QCD is not presently feasible. For this reason 
it is reassuring that, as we explained in section 2, they are both realized in all gauge theory 
plasmas with a gravity dual in the large-iV c , strong coupling limit. In this section we will 
assume that the two properties are also realized in the QGP and extract some phenomeno- 
logical consequences that might be observable in heavy ion collisions. 12 Since the heavier the 
meson the more perturbative its properties become, we expect that our conclusions are more 
likely to be applicable to the charmonium rather than to the bottomonium sector. 

An interesting feature of the energy loss by Cherenkov radiation is that, unlike other 
energy-loss mechanisms, it is largely independent of the details of the quark excited state, 
such as the precise features of the gluon cloud around the quark, etc. In the gravity description 
these details would be encoded in the precise profile of the entire string, but the Cherenkov 
emission only depends on the trajectory of the string endpoint. This leads to a dramatic 
simplification which, with the further approximation of rectilinear uniform motion, reduces 
the parameters controlling the energy loss to two simple ones: the string endpoint velocity v 
and its radial position tq. The former is just the velocity of the quark in the gauge theory, 
whereas the second roughly measures the size of the gluon cloud that dresses the quark [15]. 
In order to obtain a ballpark estimate of the magnitude of the energy loss, we will assume 
that in a typical collision quarks are produced with order-one values of tq. Under these 
circumstances the energy loss is of order unity in units of (2irT) 2 /N c , which for a temperature 
range of T = 200 - 400 MeV and N c = 3 leads to dE/dx « 2 - 8 GeV/fm. This is is of the 
same order of magnitude as other mechanisms of energy loss in the plasma; for example, the 
BDMPS radiative energy loss dE/dx = a s CpqL/1 yields values of dE/dx = 7 — 40 GeV/fm 
for q = 1 — 5 GeV 2 /fm, a s = 0.3 and L « 6 fm. Since our gravity calculation is strictly valid 
only in the infinite-quark energy limit (because of the linear trajectory approximation), we 
expect that our estimate is more likely to be applicable to highly energetic quarks at LHC 
rather than to those at RHIC. 

Even if in the QGP the magnitude of Cherenkov energy loss turns out to be subdominant 
with respect to other mechanisms, its velocity dependence and its geometric features may still 
make it identifiable. Indeed, Cherenkov energy loss would only occur for quarks moving at 

10 In some models light mesons also remain bound above T c by Coulomb-like forces [25]. 

11 An alternative possibility would be that meson states with q above some upper bound cease to exist. In 
any event, our conclusions rely on meson states existing (with a sufficiently narrow width) only up to some 
moderate q for which q > uj(q). 

12 Implications for photon production have been discussed in [28], and for deep inelastic scattering in [29]. 
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velocities v > v lim , with v lim the limiting velocity of the corresponding meson in the plasma. 
The presence of such a velocity threshold is the defining characteristic of Cherenkov energy 
loss. The precise velocity at which the mechanism starts to operate may actually be higher 
than -u lim in some cases, since the additional requirement that the energy of the quark be equal 
or larger than the in- medium mass of the meson must also be met. A related conclusion of 
our calculation in the M = 4 model is that the energy loss decreases as the velocity of the 
quark approaches the speed of light. This is in fact a universal feature of all plasmas with 
a gravity dual. The reason is that the quark velocity can approach unity only in the limit 
ro — > oo, in which the effective couplings (5.17), (5.34) and (5.55) between the quark and 
the mesons vanish because the meson radial wave functions are normalizable. In fact, in the 
case of scalar mesons the energy loss ceases completely at some subluminal velocity at which 
(5.49) vanishes. If these properties also hold in QCD then Cherenkov energy loss may be 
identifiable because it only operates in a limited range of quark velocities. 

Cherenkov mesons would be radiated at a characteristic angle cos 9 C = v lim /v with respect 
to the emitting quark, where v is the velocity of the quark. Taking the gravity result as 
guidance, v lim could be as low as v lilD = 0.35 at the meson dissociation temperature [10], 
corresponding to an angle as large as 9 C ~ 1.21 rad. This would result in an excess of 
heavy mesons associated to high-energy quarks passing through the plasma. Our estimate 
of the energy loss suggests that the number of emitted J/Y>'s, for example, could range from 
one to three per fm. This emission pattern is similar to the emission of sound waves by 
an energetic parton [30] in that both effects lead to a non-trivial angular structure. One 
important difference, however, is that the radiated heavy mesons would not thermalize and 
hence would not be part of a hydrodynamic shock wave. As in the Mach cone case, the 
meson emission pattern could be reflected in azimuthal dihadron correlations triggered by a 
high-p T hadron. Due to surface bias, the energetic parton in the triggered direction is hardly 
modified, while the one propagating in the opposite direction moves through a significant 
amount of medium, emitting heavy mesons. Thus, under the above assumptions, the dihadron 
distribution with an associated J/ip would have a ring-like structure peaked at an angle 
9 « vr - 9 C . 

A final observation is that Cherenkov energy loss also has a non-trivial temperature 
dependence, since it requires that there are meson-like states in the plasma, and therefore it 
does not take place at temperatures above the meson dissociation temperature. Similarly, it is 
reasonable to assume that it does not occur at temperatures below T c , since in this case we do 
not expect the meson dispersion relation to become spacelike. 13 Under these circumstances, 
the Cherenkov mechanism is only effective over a limited range of temperatures T c < T < T disa 
which, if T diss > 1.2T C as in [31], is a narrow interval. As was pointed out in [32], a mechanism 
of energy loss which is confined to a narrow range of temperatures in the vicinity of T c 
concentrates the emission of energetic probes to a very narrow layer on the collision geometry 
and is able to explain v 2-data at high pr at RHIC [33] . Provided that the meson dissociation 

13 This assumption is certainly correct for plasmas with a gravity dual, since the corresponding geometry 
does not include a black hole horizon if T < T c . 
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temperature T diss is not much larger than T c , the radiation of Cherenkov mesons is one such 
mechanism. 

7. Discussion 

Cherenkov emission of mesons [34] and gluons [35] in QCD has been considered before. Al- 
though some of the underlying physics is similar, the mechanism we have discussed is different 
in several respects. First, it operates in the QGP, as opposed to in a hadronic medium as in 
[34], and the radiated particles are colourless mesons, as opposed to gluons as in [35]. Second, 
the gauge/string duality provides a large class of completely explicit examples (although none 
of them includes QCD) in which this mechanism is realized and in which the energy loss can 
be calculated without further model assumptions. 

Figs. 11, 12, 13 and 14 show the rate of energy deposition into vector and scalar modes on 
the brane. All these figures share the property that the energy loss diverges as 1/Vq (as shown 
analytically in appendix D) in the limit tq — > 0. However, this mathematical divergence is 
removed by physical effects we have not taken into account. For example, for sufficiently 
large q the radial profile of the mesons becomes of order the string length and stringy effects 
become important [17]. Also, mesons acquire widths r oc q 2 at large q [36] and can no 
longer be treated as well defined quasiparticles. Finally, the approximation of a constant-?;, 
constant-ro trajectory ceases to be valid whenever the energy loss rate becomes large. 




v 



Figure 15: Energy loss into the transverse vector mode -4(„ =0 ) f° r an embedding with m = 
1.32, Rq = 1.2 (left-hand side curves) and m = 2.0, Rq = 2.0 (right-hand side curves). The con- 
tinuous curves correspond (from top to bottom) to r$ = 0.58,0.70,0.83,0.98, 1.18 (left-hand side) and 
to ro = 1.50,1.69,1.91,2.18,2.52 (right-hand side). The dotted curve is defined by the endpoints of 
the constant-ro curves. 

Figs. 11, 12, 13 and 14 also illustrate the simple dependence of the energy loss on the ratio 
m oc M q /T. Increasing m decreases the redshift at the bottom of the branes, and therefore 
increases the limiting velocity of mesons, v lim , at which quark energy loss via Cherenkov 
emission starts to operate. This means that the energy loss becomes concentrated on a 
narrower range of velocities, closer to unity, as m is increased, but the structure of the curves 
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Figure 16: Energy loss into the scalar mode ip) n _J for an embedding with m = 1.32, i?o = 1.2 
(left-hand side curves) and m — 2.0, Rq — 2.0 (right-hand side curves). The continuous curves 
correspond (from top to bottom) to r = 0.64,0.74,0.86,1.00,1.18 (left-hand side) and to r = 
1.50, 1.69, 1.91, 2.18, 2.52 (right-hand side). The dotted curve is defined by the maxima of the constant- 
r curves. 

is roughly the same up to a rescaling. This can be seen in the figures above by comparing the 
energy loss for m = 1.2 (left-hand side) and m = 2.0 (right-hand side). The concentration 
of energy loss on a narrower velocity interval is also illustrated in figs. 15 and 16, where the 
result for both values of m is shown simultaneously on the same plot (for the transverse vector 
and scalar modes). 

In this paper we have concentrated on the case N { = 1, i.e. we have assumed the presence 
of a single heavy flavour. Consider now a theory with multiple heavy quarks, such as QCD 
with c and b quarks, for example. In the string description this corresponds to a situation with 
iVf > 1 D-branes. If all quarks have identical masses (and R-symmetry quantum numbers) 
then the D-branes are all coincident and their worldvolume theory is described by a non- 
Abelian U(N f ) theory, corresponding to the fact that mesons come in multiplets that 
transform in the adjoint representation of U(N t ). A quark with flavour / then may emit any 
of the iV t mesons with flavour //, with / = 1, . . . , N { . Under these circumstances the energy 
loss is enhanced by a power of iV f . 14 Note that if the resulting meson has f f then the 
quark must change flavour / — >■ f in the emission process, as shown in fig. 17(left). A process 
in which the quark does not change flavour is also possible, as shown in fig. 17(right), but 
this requires the emission of at least two mesons and is therefore further suppressed at large 
iV c . 

Consider now the opposite situation, more analogous to QCD, in which the heavy quarks 
have different masses, so that the D-branes in the string description do not overlap. This is 
depicted in figs. 18, 19 and 20, which show the string description of the emission processes of 
figs. 17(left) and 17(right), respectively. Figs. 18 and 19 correspond to emissions with / ^ /, 
whereas fig. 20 describes the emission of an //-meson. In this geometric picture the necessity 

14 This corrects the corresponding statement in ref. [1]. 
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Figure 17: Meson emissions by a quark in a theory with multiple flavours. 



of a two-meson emission in the case / 7^ / in order to preserve the quark flavour is due to the 
fact that the string must break twice in order to stay attached to the same brane - see fig. 19. 
Since string breaking is suppressed at large iV c , this process is subleading with respect to one- 
meson emission. In any case, since an open string must always have its endpoints attached to 
a brane, the emission of one or multiple mesons with / 7^ / by a quark of flavour / requires a 
tunneling process in which the string fluctuates and touches the /-brane, as shown in figs. 18 
and 19. The amplitude for this process can be studied semiclassically provided the distance 
between the two branes is sufficiently large compared to £ s , but it is far from straightforward 
to calculate [37]. In addition, in the present context the calculation would require a precise 
specification of the string profile. On general grounds, however, one may expect the amplitude 
to be exponentially suppressed, since it requires a large string fluctuation that is classically 
forbidden. Note that the same exponential suppression applies to the emission of a large //- 
meson by an /-quark, as shown in fig. 20. This is the reason why we neglected this process 
in our calculation of energy loss. More precisely, our calculation can be seen as accounting 
for this process in the limit in which the size of the emitted meson is so small that it requires 
quantization of the resulting string. In this limit there is no exponential suppression, and the 
emitted string must be described as a field propagating on the brane. 

We close with a comment on a possible extension of our work. In this paper we have 
focused on the energy loss of quarks attached to branes that sit outside the horizon in a 
Minkowski embedding. It would be interesting to study the energy loss in the case of black hole 
embeddings, which describe light quarks. In this case no stable quark-antiquark bound states 
exist in the plasma, which in the string description corresponds to the fact that excitations 
on the brane are characterized by quasinormal modes with complex frequencies. It would 
be interesting to explore whether the emission of quasinormal modes could also lead to a 
significant energy loss. 
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Figure 18: String description of the emission of an //-meson with / 7^ /. The quark changes flavour 
in the process. 




Figure 19: String description of the emission of two mesons with / ^ /. The quark does not change 
flavour in the process. 




Figure 20: String description of the emission of an //-meson by an /-quark. The quark does not 
change flavour in the process. 
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A. Boundary conditions at the string endpoint 

The action for the string may be written as 

S = -T Btring J drda^g l - g a ^d a X M 8^X N G M N • (A. 1) 

In this form of the action g is an independent worldsheet metric and X m (t, a) specify the 
embedding in spacetime of the string worldsheet. Since g appears undifferentiated, it can be 
eliminated from the action through its equation of motion. This equation implies that 

g a p = d a X M d p X N G M N , (A.2) 

i.e. that g is the induced metric on the worldsheet. Substituting this into the action one 
obtains the familiar Nambu-Goto action 

S = _T^ j dTda ^- g (A 3) 

An alternative way to proceed, which is more convenient to elucidate the boundary 
conditions at the string endpoint, is to choose the so-called conformal gauge. This means that 
one uses the reparametrization invariance of the string action to ensure that the worldsheet 
metric is conformally flat, i.e that g a p = Q 2 (t, <j)r) a p- (In addition, Weyl invariance may be 
used to ensure that Q, = 1.) In this gauge the action becomes 

S = -T string J drda 1 - ri a Pd a X M dpX N G MN . (A.4) 

Variation of this action with respect to the embedding coordinates yields a bulk term pro- 
portional to the equation of motion, 7] a ^d a dpX M = 0, plus the boundary term 

(A.5) 

bdry 

which is integrated over the string boundary. The equation of motion must be supplemented 
by the constraints associated to the gauge fixing of g, which take the form 

A. ■ A. = Lrj^jyvV — U , 

X 2 + X' 2 = G M n(x m X n + X M 'X N '} =0. (A.6) 

The boundary conditions follow from the requirement that the boundary term vanish. This 
may be achieved by imposing either a Neumann boundary condition, X M/ | bdry = 0, or a 
Dirichlet boundary condition, 5X M \ hdiy = 0. If all coordinates satisfy Neumann boundary 
conditions, then the second constraint immediately implies that AT 2 | bdry = 0, namely the 
familiar condition that the endpoint moves at the speed of light. Suppose however that the 
string is attached to a Dp-brane, and let the first p+1 coordinates X a be coordinates along 
the brane directions, and X A be coordinates orthogonal to the brane. Then by definition 



dr 



vM' r yi\ n 
A OA (jMN 
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we must choose Neumann boundary conditions for X a and Dirichlet boundary conditions for 
X A : 

^ a '|bdry = 0, X A \ biry = . (A. 7) 

The Neumann boundary condition on X a implies that the string ends orthogonally on the 
brane, since the vector tangent to the string at its endpoint, X M '\ hdry = 0, has no components 
along the brane. Substituting both boundary conditions on the second constraint equation 
one finds that ^ 

X |bdry = ~ X | bdry ^ 0, (A-8) 

which means that the endpoint moves along the brane at a speed lower than or equal to the 
local speed of light. 

B. Energy loss formula 

A charge can be defined for each of the isometries of the brane's worldvolume metric (3.6). For 
concreteness, in this section we will focus on the four-momentum associated to translations 
in the gauge theory, which is generated by the set of Killing vector fields k^ = d^. 

Let £ be a spacelike 7-surface in the brane's worldvolume, which we take to be a t = const, 
surface, and n = dt/^—goo the future-pointing unit normal to S. The momenta are then 
given by 

P» = j d 7 x^gT p n a T ab fcj = J d 7 x ^g-fi a T afl , (B.l) 

where g sp is the spatial metric on S. Since the time-space off-diagonal components of g vanish 
we have that 

g 00 = 1/500 , = ^m^gTv ■ (B.2) 

Using these relations may be rewritten as 

P» = f d 7 x^gV p —^T^ = - [ d 7 x^gT° (B.3) 
Jt V~9oo Jt 

For fj, = 0, these formulas give the energy on the brane: 

E = P = [ d 7 x^ p —^T 00 = - [ d 7 x^T° >0. (B.4) 
Jt V~9oo Jt 

Note that this is non-negative because Too > 0. 15 

Consider now the brane's worldvolume V as shown in fig. 21. Si and X2 are spacelike 
hypersurfaces at times t\ and £2, respectively, and £00 is a timelike hypersurface at spatial 
infinity. Applying Stoke's theorem we then have 

I V~~~g~V a T a ^ = f Jg- n a {l) T a ^+ ! ^ n\ 2) T a(1 + f ^ n a {oo) T aix . (B.5) 

15 For example, in flat space eq. (5.6) gives T 00 = E 2 /2 + B 2 /2 > 0, where E t = F i0 , E % = F%, E 2 = E i E i , 
and 2B 2 = FaF ij . 
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The unit normals must be taken inward-pointing if they are time- 
like, and outward-pointing if they are space-like. Therefore we 



have 
yields 



-n 



(2) 



n and nf % 



n, so eq. (B.5) with fj, = 



-gV a T c 



oO ■ 



(B.6) 



where E\ 2 is the energy contained in Si 2 and F^ = — L ^/g^n a T a Q 
is the energy flux that has escaped through between t\ and 
£2- 16 The left-hand side of eq. (B.6) is the total energy deposited 
on the brane by the source. Since f v = J dtd 7 x, if t\ and £2 are 
innitesimally close we obtain 

dE _ 
~dt ~ ' 




d 7 x^ 



-gV a T c 



aO ■ 



Figure 21: Brane's world- 
volume, to which Stoke's the- 
orem is applied. 

(B.7) 



C. High-momentum radial wave functions 

As shown (on the left-hand side of) figs. 6, 8, and 9, the high-momentum radial profiles of 
the different modes are concentrated near the tip of the brane (r = 0). In this region of high 
q and small r it is possible to find analytic expressions for the radial profiles [17]. In this 
appendix we shortly review this computation and extend it to vector modes. 
Following [17] we introduce the coordinate z which fulfills 



dz 
dr 



9r 



'2/(1 + R 2 ) 



(C.l) 



-500 V p 4 f 2 

In terms of this new coordinate, the different modes *$> a = {fa, fa, A, 3>} satisfy a differential 
equation of the generic form 



d z [a a (z)d z * a ]+a a (z) 



f 2 

u? - Lq 2 
P . 



b a (z) (m a )H Q 



0. 



(C.2) 



where m a = {mix, 77122, 0,0} and the different coefficient functions are given by 



- 9 9 00 , 

-gg 00 g n 



-AiiSooA 
l^g- 



(C.3) 



Via the simple rescaling ip a 
form 



'aP- ty a , the set of eqs. (C.2) can be written in the Schrodinger 



-d'i V> Q + V a ^ a = co 2 ip a , 



(C.4) 



1 The minus sign in the definition of the flux comes from the fact that the energy current in the a-direction 
is given by — T a o, as can be seen from the continuity equation. For example, in flat space in the absence of 
sources the continuity equation d a T a o — yields doToo — diTm = 0. 
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where the potential is given by 



V a (z,q)=q 2 ^+ z ' + — (m a y. (C.5) 

This potential is a complicated function of z which is only known numerically because it 
depends non-trivially on the brane embedding. However, since the radial profiles are concen- 
trated near the tip of the brane, we focus on the small- z expansion of the potential (C.5). In 
the limit z — > 0, eq. (C.l) reduces to 17 

Z= Rt-l r ' (C - 6) 

and the small- z potential is given by 

3 \ 1 2 2, 1 2 „2 n 2 



V a (z, ?) = + n** + 2) J ^ + + ^ W , (C7) 

where i a = {0, 0, 0, 1}, v lim is the meson limiting velocity (3.11), and 

2 lQRl{Ri-lf{l + Rl) 

n = OTW • (C ' 8) 

As claimed, this potential has a minimum at z oc 1/y/q which means that, at least for the 
lowest modes, the wave functions are concentrated at small z. Different meson excitations 
correspond to different states of the four-dimensional harmonic oscillator (C.7). The eigen- 
functions and eigenvalues are given by 

= vf iin q 2 + qn(2n + 2 + e a ), (C.9) 
V£ = M a zl +r L i i a+l ) Q^ 2 ) exp , (CIO) 

where is the generalized Laguerre polynomial and N a are normalization constants de- 
termined by the requirements 



/ 



dz^(z)^(z) = S mn . (C.ll) 



These normalization conditions coincide with those in eqs. (4.9), (4.20) and (4.32). 

For future use, we provide explicit expressions for the lowest excitations, which correspond 
to n = 0. The normalization constants are given by 



2 



17 To derive this expression we have used the small-r expansion of the embedding, R(r) — Ro + 
r 2 /(R (R 8 -l)). 
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and the radial wave functions take the form 



where f3 a = lim r ^o Va"'- 



o3/2 o3 

P' = P 2 = : 4 V (C14) 
(i? 4 -l) 3/2 

, 2^(1 + ^1/2 



/3 4 = \° =^7779 - (C16) 



4 n\3/2 



^-1 



2fl§ (l + jR 4) 



4\l/2 ' 



D. Energy loss at small r 



As the quark position approaches the tip of the branes, ro — > + , the maximum velocity of 
the quark approaches the meson limiting velocity t> lim from above. As a consequence, the 
quark and meson dispersion relations cross at go S> T. In fig. 2 this means that the dotted 
vertical lines move to the right. For a fixed ro the smallest value of the crossing point, qy in , 
is attained at the maximal velocity of the quark, w max . Using the high- momentum dispersion 
relation (C.9), the crossing momentum at an arbitrary quark velocity v is determined by the 
condition 

fiim<7o + 1 + n + — =vq , (D.l) 

V 2 J u lim 

which leads to 

qo=(l+n+ e pi—^—. (D.2) 



The maximal velocity for a quark at ro is 

^maxOo) = 



goo 

9u 



tt 2 

~ flim + Zq , (D.3) 

r 



where we have expanded to leading order in ro and used the definition (C.6) of z. Substituting 
this value of v in eq. (D.2) we find that the minimum crossing point is 

*-(»+»+T)i5r (D ' 4) 

Thus, the energy loss can be reliably computed with the approximate solutions of appendix C. 
Since the energy-loss formulas are different for each mode, we will address them separately. 
Furthermore, we will focus on the lowest state of each mode. 
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D.l Scalar mesons 



Since in the high-momentum, small-r region both scalar modes have the same radial profile, 
the energy loss into scalar mesons in this limit will also be the same. Expanding the energy 
loss formula (5.57) and using the radial profiles (C.13) we obtain 



dt 2Q 3 
After integration, this yields 



Rq 



2Rl 



dq 
2^ 



<f exp ( 



(D.5) 



dE. 



e 2 (l + ^) 2 (^) 



scalar 



^ 2 r v 2_ 



dt 



4ttO, 



2RI 



X z o) ~ v 



21 



48 + 24Qq$ in z% + 6tt 2 (qfi™) 2 z$ + 3 (gg lin ) 3 ^ 



exp 



. (D.6) 



As expected on general grounds, and in agreement with figs. 13 and 14, we see that the energy 
loss vanishes both for v — > v max and for v — > i> lim . The former is due to the fact that the factor 
(5.49) vanishes in this limit. The latter is implemented by the fact that c/JJ" n diverges as 
v — > fn m , which in turn is a manifestation of the Heaviside theta function in eq. (5.57). As 
reflected in figs. 13 and 14, the maximum energy loss occurs at some intermediate velocity such 
that u lim < v int < f max . Although v int is not easy to compute, we know that t> max — v lim oc z 2 . 
It then follows that also v max — v int oc Zq and therefore that the maximum energy loss diverges 
as I/zq. 

D.2 Transverse vector mesons 

The energy lost into these modes is given by eq. (5.23) which, utilizing (C.13), leads to 



dE t , 



dt 



e 2 v 



r 



dq 3 
2^ 6XP 



1 



QqzQ 



1 



v 



2 

lim 

2~ 



2 A) 

qv z J 



(D.7) 



Unlike for scalar meson emission, it is easy to see that in this case the energy loss is a 
monotonically growing function of the velocity. Thus, the maximum value of the energy loss 
is attained for v = v max and, to leading order in zq, it is given by 



dE u 



2 76 exp(— 4) 



VlimZQ 



dt AttQ 3 

This shows the same divergence for small ro as in the case of scalar modes. 



(D.8) 



D.3 Longitudinal vector mesons 

Inserting (C.13) into the expression for energy loss (5.40), and to leading order in zq we obtain 



dE, 



long 



e 2 2ft 3 1 



dt 



W 4 



2^ ' ' 



VLqzQ I exp 



QqzQ 



(D.9) 
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As in the transverse vector case, the maximum energy loss is attained at v = f max . Using 
(D.4) and setting £ a = 1 (see the definition below eq. (C.7)) the maximum energy loss is 



dE long _ e 2 1 632exp(-6) 
dt ~4^h~(f^f v, im 4 



(D.10) 



Again, this diverges with the same power of zq as in the cases of transverse and scalar modes. 
Note also that since (3 A = 1//3 3 and the numerical factors in eqs. (D.10) and (D.8) are similar, 
we find that in this limit the energy lost into transverse and longitudinal vector modes is 
comparable. 

Finally, we may comment on the mass dependence of the energy lost into vector mesons. 
The maximal energy loss is proportional to 



dE 
lit 



Rt{Rt-lf 1 



oc " v " > 4. (D.ll) 
(Rt + l) rl 



The fact that this is a growing function of Rq means that the energy loss increases if the 
quark mass is increased while keeping all other parameters such as vq, v, etc. fixed. However, 
this should not be necessarily taken as an indication that Cherenkov energy loss in a real HIC 
experiment increases as the quark mass increases, because the 'preferred' or 'mean' values of 
these parameters with which a quark is produced may themselves depend on the quark mass. 

E. Low-temperature limit 

As T — > with all other scales fixed, the redshift at the bottom of the branes decreases, 
and so the limiting velocity of mesons, v Hm , approaches unity. In turn, this means that the 
momentum qo(T) at which Cherenkov radiation turns on diverges as T — > 0, much in the 
same way as in the case ro — > studied above. The purpose of this section is to estimate 
qo(T) in the low-temperature limit. As we will see, the product qo(T)T 2 remains finite as 
T — > 0, which makes this limit harder to study than the ro — > limit. 

As T — > with fixed quark mass, the parameter m oc MJT controlling the asymptotic 
position of the branes, R(r — > oo) = m, diverges. For this reason it is convenient to introduce 
rescaled coordinates 

x» = — , {f,R,p} = V~e{r,R, P }, (E.l) 
v e 

where e = 1/m 2 and we are interested in the limit e — > 0. The rescaling of the x-coordinates 
is chosen so that the induced metric on the D7-brane takes the same form as that in (3.6), 
i.e. ds 2 = L 2 ds 2 (g) with 



.2 



ds 2 (g) = P - 

where now 



f 2 

—dt 2 + fdx 2 
. f 



?2\ -2 



+ { ± ± ^- 1 df 2 + - 2 dn 2 , (e.2) 



f(p) = l-y 4 , /(P) = l + ^. (E.3) 
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In this new set of coordinates the horizon is located at p hor = y/e, whereas R(r — > oo) = 1. In 
fact, in the limit e — > the gravitational pull of the black hole becomes very small and the 
brane bends very little. Inserting the ansatz R(r) = 1 + 5R(f) in eq. (3.7) and linearizing in 
5R one finds that SR = 0(e 4 ). Since we will work to order e 2 we will neglect this correction. 
Furthermore, it is easy to see that the equations of motion for the transverse (4.6), longitudinal 
(4.18) and scalar (4.30) modes are unmodified by the rescaling (E.l) provided we rescale the 
momentum in a consistent way, i.e. {(2>,q} = \fe{uj, q}, so that we leave the product (ujt — q-x) 
invariant. For ease of notation, in the following we will omit the hat symbol. 

Following appendix C, eqs. (4.6), (4.18) and (4.30) can be rewritten in the Schrodinger 
form (C.4). The potential (C.5) is in general a complicated function, but it simplifies in the 
limit e — > 0. Indeed, in this case eq. (C.l) can be integrated explicitly with the result 

3 / r(5 + 3r 2 ) \ 

z(r) = \/2arctanr+ — = 3arctanr+^ / e 2 + C(e 4 ) , (E.4) 

8^ V (1+r 2 ) 2 J 

so that z max = z(r — > oo) ~ ir/y/2. The potential can then be written as 

V a (z,q) = V 1 (z)q 2 + V%(z), (E.5) 

where 

V^z) = l-4cos 4 (-^)e 2 + 5 (z)e 4 + 0(e 6 ), (E.6) 
2sm z h/2z l + cos(v22) 



with A a = {0, 0, 1, 1}, B a = {-6, -6, 0, 8} and C a = {9, 9, 2, 6}. The functions g(z) and h(z) 
are smooth and bound, and their explicit form will not be needed. 

Before proceeding further let us clarify one point. In the limit z — > z max , V\{z) shows a 
negative divergence. This is related to the fact that the radial profile ^ 4 (r) = <E>(r) does not 
vanish near the boundary. This may seem counterintuitive, since generically one expects a 
'confining' potential near the boundary (the AdS 'box') which is partly responsible for the 
discreteness of the spectrum. However, note that the physical electric field E is related to $ 
through eq. (4.16), which may be written as E = a 4 (z)d z &. In the limit r — > oo, the factor 
f 2 /f 2 in eq. (C.2) approaches 1, so one can differentiate this equation with a = 4 and use the 
fact that a A (z) = l/a 3 (z) to show that the asymptotic form of the equation for E is identical 
to that for the transverse mode \I /3 = A. 

Let us now return to the potential (E.5). As explained above, we are interested in the 
limit e — > and q > qo{t) — > oo. In particular, we wish to determine whether in this limit 
the product eqo goes to zero, remains finite, or diverges. We will establish that the product 
remains finite by showing that the other two possibilities lead to a contradiction. 

Consider first the possibility that ego — > oo as e — > 0. Then the potential takes the form 

2 A a cos 2 (^)+B a cos( v / 2z) + C Q (m°) 2 4 /z \ 22 

V (z, q)-q 2 ~ ' q / i — \ + 7-p-r - 4 cos -p e 1 . 

2sin 2 ( v / 2^) l + cos( v / 2^) \V2J 

(E.8) 
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where the last term dominates everywhere except near the endpoints z = and z max , at 
which the order-one part of the potential, given by the first two terms on the right-hand side, 
diverges. As an illustration, fig. 22 shows the potential for the a = 3 - mode for several values 
of eq. We see that the potential develops a minimum at a small value of z as eq becomes large. 
This allows us to consider an expansion for 
small z in order to obtain the energy levels 
associated to the potential. In this way we 
reproduce the result (C.9), approximated for 
e — > 0. For example, in the case a = 3 and 
n = the dispersion relation takes the form: 
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(E.9) 

The value of qo is determined by the condi- 
tion co(qo) = vqo, where v is the quark velocity, 
which depends on the quark position through 




0.0 



Figure 22: Potential for a — 3, in the limit 
e — > 0, q — > oo and eq fixed. From top to bottom 
the curves correspond to eq = 4, 8, 12, 16, 20. 
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The result is 



-10 + VE 



20 + r2(2 + r§) (40 + 21r2(2 + r 2 ,)) 

(1 + 7-2)2 



-1/2 



(E.ll) 



This expression is not parametrically large for any value of ro except in the limit ro — > 0, which 
was considered above and is unrelated to the zero-temperature limit under consideration here. 
In particular, for small ro eq. (E.ll) yields eqo = l/ro 2 + 3/2 + O(t-q), whereas for large ro 
one finds eqo = 2.01 + 1.99/Vg + ©(I/Vq). We therefore conclude that eqo(e) remains finite as 
e — > 0, in contradiction with our initial assumption that eqo(e) — > oo in this limit. 

Consider now the opposite possibility, i.e. that eqo — > as e — > 0. In this case the last term 
in the potential (E.8) is small and the energy levels can be determined using perturbation 
theory with the result 



+ A n - w n e 2 q 2 + 0(eV 



,4„4n 



(E.12) 



Here, A n are the eigenvalues of the problem in the absence of the perturbation, and —uj n is 
the (negative) first-order correction given by the expectation value of the perturbation in the 
n-th eigenstate, (n\ — 4cos A (z/\/2)\n) . The key point is that \ n and w n are independent 
of e and q. In order to find the crossing point, we need to solve w 2 ((/o) = v 2 q$. Since the 
minimum possible value of qo corresponds to the maximum quark velocity, v = 1, we find 




(E.13) 
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which is in contradiction with our initial assumption that eqo — > 0. 
We therefore conclude that the combination 



1 q ( V\T \ q /2vrT\ 3 q _ ^ 



remains finite in the limit T — > 0, where we have reinstated the hat and we recall that q is 
the physical, dimensionful momentum. We see that in this limit q(T) T 2 ~ M q 3 /A 3 / 2 ~ M 3 ea , 
which remains finite in the low-temperature limit, as we anticipated. This means that the 
relevant potential in this limit is (E.8), where all terms are of the same order. This makes the 
problem harder than that associated to the limit ro — > 0, and we have been unable to find 
analytic expressions for the corresponding eigenfunctions and eigenvalues. 
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